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SUMMARY 

A general steady through-flow theory of nonviscous fluid in turbo- 
machines,' 1 of arbitrary hub- and casing -wall shapes with subsonic or 
supersonic velocity is presented. The theory is applicable to both 
direct and inverse problems and is derived primarily for use in turbo- 
machines having thin blades of high solidity with a simple approximate 
correction factor for blade-thickness effect. Through the use of the 
stream function, the continuity equation and the equation of motion in 
the radial direction are combined to form a principal equation for the 
present problem. The principal equation contains some terms that are 
either prescribed or to be determined by other equations defining the 
problem. Two forms of the principal equation are obtained for the two 
main groups of current compressor and turbine design in which the vari- 
ation of tangential velocity and the variation of the ratio of tangen- 
tial to axial velocity throughout the blade region are given. When the 
tangential velocity is given, the principal equation is elliptic or 
hyperbolic, depending on whether the meridional velocity is subsonic or 
supersonic, respectively. When a relation between the tangential and 
the axial velocity is given, however, the principal equation becomes 
hyperbolic when the relative velocity is supersonic. A general method 
of solution for both the elliptic and the hyperbolic case is outlined. 
Specific applications of the theory to several common types of com- 
pressor and turbine employing free-vortex, symmetrical-velocity-diagram 
solid-rotation- type, nontwist ed-blade, and radial -blade -element designs 
are discussed. 


INTRODUCTION 

tfith the increasing use of velocity diagrams other than free- 
vortex type, low inlet hub -tip -radius ratios, and high velocity of flow 
the problem of three-dimensional flow in axial-flow turbomachines 
becomes more and more important. This problem is treated by Traupel, 
Meyer, and Marble (references 1 to 3) for incompressible fluid. 
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Compressible flow is treated in reference 4, in which methods for 
limiting solutions of zero and inf inite -hlade aspect ratio are obtained 
and a step-by-step method, as well as a simpler method based on an 
approximate knowledge of the shape of streamlines, is given for finite- 
blade aspect ratio. In reference 5, Reisner gives a method of blade 
design for compressible flow with the shape of neither hub nor casing 
wall specified in advance. The problem of supersonic flow in impellers 
of a given casing and blade shape is currently being investigated by 
Arthur W. Goldstein of the NACA Lewis laboratory. 

The analysis made at the NACA Lewis laboratory and presented 
herein proposes a unified theory that is applicable to both direct, and 
inverse problems for both subsonic and supersonic flows in compressors 
and turbines of arbitrary hub- and casing-wall shapes. 

Equations of motion and energy for unsteady three-dimensional flow 
of a nonviscous compressible fluid are expressed in terms of some con- 
venient quantities for analyzing flow in turbomachines. Entropy changes 
due to heat transfer in a cooled turbine and due to shock wave in 
supersonic flow can easily be accommodated. The condition under which 
irrotational-flow analysis is correct is also discussed. 

The general equations are then simplified for steady through-flow 
in turbomachines having thin blades of high solidity. It is shown 
that, in the direct problem, just enough equations exist to determine 
all the variables; whereas in the inverse problem, after the inclusion 
of the integrability condition for the blade surface, either one vari- 
able or a relation between several variables can be specified by the 
designer. 

In the solution of the problem, the continuity equation and the 
equation of motion in the radial direction are combined into a principal 
equation through the use of the stream function. This equation 
involves some terms that are either given or to be determined by other 
equations defining the present problem. Two forms of the principal 
equation are obtained for two main groups of designs in which the vari- 
ation of tangential velocity and the variation of the ratio of tangen- 
tial and axial velocities are prescribed by the designer. The crite- 
rions of whether the principal equation is elliptic or hyperbolic are 
obtained for both groups. 

A general method of solving the set of equations for both the 
direct and the inverse problems is then described for turbomachines of 
arbitrary hub- and casing-wall shapes and for either an elliptic or 
hyperbolic principal equation. 
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The dimensionless forms of the principal and other equations for 
some typical designs are given. A simple approximate correction for 
blade thickness is also given. 


SYMBOLS 


The following symbols are used in this report: 


a 


velocity of sound 


m i differentiation coefficients in equation (70) used to multi 

n J ply function value at point Xj to give m LLL derivative 

at x^ using polynomial of n^ 1 degree 


b,f,J,K, functions of r and z 
L,M,N 


C 


constant 


_D_ 

Dt 

F 

G 

H 

h 

I 

k 

l 

m 


specific heat of gas at constant pressure and volume, 
respectively 

differentiation with respect to time following motion of 
fluid particle 

mth derivative of q 

blade force per unit mass of fluid 

Green's function 

total enthalpy per unit mass of fluid, h 
enthalpy per unit mass of fluid 

h + T " T = H " v u r) 

thermal conductivity 
characteristic curve 
order of derivative 
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n 

P 

Q 


q. 

R 

m i 
n* 

r, 6,z 
r,q>,z 
s 

T 

t 

U 

u 


V 

W 

x 




r 

r 


degree of polynomial 
pressure 

heat added to fluid particle along path of motion per unit 
mass per unit time 

dependent variable 

gas constant 

remainder term of m^h derivative at point xp obtained by 
using polynomial of n™ 1 degree 

cylindrical coordinates relative to stator 

cyl in drical coordinates relative to rotor 

entropy per unit mass 

temperature 

time 

velocity of blade at radius r 

internal energy per unit mass with 0° absolute as base 
temperature 

absolute fluid velocity 
relative fluid velocity 
independent variable 
W u 

arc tan rp- 
w z 

grid spacing in r- and z-directions, respectively 

ratio of specific heats 
average value of y 
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independent variables 

dr 

slope of characteristic curve* ^ 

W u 

tan p = 

value of x between x Q and x n 
mass density- 
stream function 
angular velocity of blade 

Subscripts : 

b trailing edge of guide vane 

c leading edge of rotor 

d trailing edge of rotor 

e leading edge of stator 

h hub 

i inlet 

j any station 

m meridional 

o refers to position where blade element is radial or F r = 0 

r,u,z radial, circumferential, and axial components 

T total state 

t at tip or casing 

Superscripts : 

a,b, grid points 

* dimensionless values 


X 

V 

P 

t 

0) 
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GENERAL BASIC EQUATIONS 

The three-dimensional compressible flow of a nonviscous gas 
through a turbomachine is governed by the following set of basic laws 
of aerothermodynamics: From the principle of conservation of matter, 

the equation of continuity is 

( 1 ) 


+ V-(pw) - o 


or 


V •¥ + 


D(log e p) 
Dt 


(la) 


For a blade rotating at a constant angular speed Newton' s second 

law of motion gives 


DW 

Dt 


co 2 r + 2c0XW = 


D V 
Dt 



( 2 ) 


The first law of thermodynamics may be written as 


v* + v pfp" 1 ) = Q 

Dt P Dt Q 


( 3 ) 


where u is related to the gas temperature T by 


Du _ DT 

Dt “ Dt 


( 4 ) 


an d Q is given by the following equation if only conduction is con- 
sidered: 


Q = p"V(k7T) 


( 5 ) 


For the range of gas temperature and pressure encountered in ordinary 
turbomachines , p, p, and T are accurately related by the following 
equation of state 


• p = RpT 


( 6 ) 
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Although the flow of gas through the turbomachine is completely 
defined by the preceding equations together with the known variations 
of c v and k with temperature and the given boundary and initial 
conditions , in reference 4 it is found more convenient to express the 
state of the gas in terms of entropy, total enthalpy, or a quantity I 
of the gas, in addition to its velocity components. These quantities 
are defined as follows: 

( 7 ) 


( 8 ) 

I = h + ^- ^ = H- aj(V u r) (9) 


T ds = du + p d 


(p- 1 ) 


H 


h + 


V 


where 


h = u + pp 


-1 


( 10 ) 


When equations (4), (6), (7), (10) , and the relation, 

R — Cp - Cy = (y-l)cy 


which follows from equations (4) , (6) , and (10) , are used, there are 
obtained 


T ds = dh 


dp 


T ds = — d| 

r-i 


(?)- 


dp 


d (f) = FI d ( loge p) “ i(log e p) 
d (D = Fl d(l ° g e T) ■ d(log e p) 

and the equation of continuity can he written as 

1 D(log e T) D 


V -W + 


T-l Dt 


Dt 


(i) = 


(ID 

(11a) 

( 12 ) 

(12a) 

(13) 
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From equations (9) and (ll) , 

- yp + i VW 2 - co 2 ? = VI - TVs 
P £ 

With this equation and the relation 

1 - g? + - I + | VM 2 - «x(vxw) 

the equation of motion (2) can he -written as 

|| - WXfrXW) + 2?XW = -VE + TVs (14) 

An alternate form of equation (14) that involves the vorticity of the 
absolute motion is obtained as follows: Usi.ng the cylindrical coor- 

dinate system with the z-axis parallel to co yields 

V = W + coxr (15) 

VXV = VX¥ + VX (cox?) 


but 


VX (coxr) = (r*V)co - (co*V)r + co(V'r) - r(V*co) = 2co 

Hence, 


VXV =VXW + 2co 


(15a) 


This relation can also be seen from the following expressions of rela- 
tive an d absolute vorticity expressed in terms of cylindrical coor- 
dinates r,cp, z and r,0,z, which refer to the rotor and stator, 
respectively: 


i dW„ dW.. 

(vxw) r = - - 'ST 


(7xw) 


dw r dw z 

u ~ ^z 5r~ 


(VXW) Z 


1 ^( w u r ) 1 ^W r 

r 5r r ctep - 



(16) 
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and the relation 


dV„ dV.. 

(VXV) r = r ST " 5T 

, _ x 5v r Bv z 
(yxv) u = ST " SF" 


S(V u r) 


(VXV) Z - i 


1 Sv r 

r d9 


a(V u r) a(¥ u r) 


"5F 


"5F 


+ 2cor 


(17) 


(18) 


When equation (15a) is used; the alternate form of equation (14) is 


|| - Wx(7XV) = -VI + TVs 


(14a) 


The use of equations (2); (9), and (ll) yields 


DI = T g? + i ^ + i(w-v)p + if-S - 5 - du 


Dt 


Dt 


Dt 


Dt 

= T i + ^st + w-(» 2 ?- fi») - o-i 


Dt p at 


(19) 


Hence, energy equation (3) can he written as 


_ rp Ds _ DI 1 dp 
Dt ” Dt ” p St 


Q = T 


( 20 ) 


The preceding equations lead to several important general considerations. 

If the blade rows are not placed too close together, the pressure 
of gas at a fixed point relative to the blade can be taken as constant 
with respect to time. Consequently, according to equation (20), the 
entropy and the quantity I of the gas stay constant along its rela- 
tive streamlines for adiabatic flow. The constancy of I means that 
the rate of change of total enthalpy along the streamline is equal to 
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the angular speed of the tlade multiplied hy the rate of change of 
angular momentum about the machine axis of the gas particle along its 
streamline, which is the well-known Euler turbine equation usually 
derived under less general conditions. In a cooled turbine where the 
heat transfer may be large, the rate of change of s and I along 
the streamline can be obtained by using equation (20) . Again, for 
steady relative flow, equation (14a) shows that the vanishing of abso- 
lute vorticity requires both gradient I and gradient s to vanish 
or the difference between VI and TVs to vanish. When both gradient 
I and gradient s are zero upstream of the blade row and the flow 
is adiabatic, s is uniform in passing through the blade rowj p is 
then a function only of p (according to equation (12)), and conse- 
quently, according to Kelvin's circulation theorem, the absolute vor- 
ticity will remain zero in passing through the blade row and the flow 
can then be treated on the basis of irrotational absolute flow. 

For flow through a stationary blade row, equation (l4a) becomes 

f - VX(V*V) = -VH + TVs (14b) 

which agrees with similar relations previously obtained by Vazsonyi 
(reference 6), and Hicks, Guenther, and Wasserman (reference 7). It 
is interesting to see that, for relative flow in a rotating blade row, 
VX(V*V) becomes WX(VXV) and H becomes I. 

When it is assumed that the fluid enters the machine with uniform 
H and s and zero vorticity, the adiabatic flow through the inlet 
guide vanes can be treated on the basis of irrotational flow. When 
the guide vanes impart a radial variation of the tangential velocity 
of the fluid downstream of the vanes similar to that in a potential 
vortex, the circulation is constant along the blade span and the fluid 
maintains a uniform H and s and a zero vorticity going into the 
following rotor blade row. If the rotor blade row is situated far 
away from the inlet guide vanes, the fluid enters the rotor with a 
uniform I in the circumferential direction as well as in the radial 
direction, and the flow through the rotor blades can again be treated 
on the basis of zero absolute vorticity. If the rotor is close to the 
guide vanes, the fluid entering the rotor blades is circumferentially 
nonuniform, which condition is balanced by the unsteady term of the rel- 
ative velocity, and the flow through the rotor blades should theoreti- 
cally be treated on the basis of unsteady flow with zero absolute 
vorticity. When the guide vanes impart a radial variation of tangen- 
tial velocity of the fluid downstream of the vanes different from that 
in a free vortex, however, the circulation varies along the span of 
the guide vanes, thereby shedding vortices from the trailing edge to 
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the fluids downstream, and the fluid enters the following rotor blades 
with a uniform s and H, a nonuniform I, and a nonzero value of 
absolute vorticity. Consequently, the flow through the rotor blade ■ 
row cannot be treated on the basis of zero absolute vorticity. 


From the preceding discussion, the choice of H or I and s as 
the basic thermodynamic variables of the gas is apparent. Compressor 
and turbine rotors are usually designed to add or subtract the same 
amount of energy radially to or from the gas; hence, H is usually 
radially constant throughout the machine. When the circumferential 
velocity pf gas upstream of the blade row is zero or varies inversely 
with radius, I is then radially constant throughout the machine when 
the heat transfer is zero or is uniform radially. These facts will be 
utilized in the following developments. 

The continuity equation (13) can also, be put into a form contain- 
ing H or I by use of a constant value of y, that is,_ y for the 
range of temperatures involved in the process. By using y, equation 
(12a) can be written as 


d(log e p) = d 





( 21 ) 


Integrating from the inlet total state yields 



where As* 
form: 


s-s 


T,i- 


The continuity equation then takes the following 



a 

St 


( 23 ) 
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The solution of the three-dimensional-flow problem thus consists 
of the simultaneous solution of equations (l) or (23), (14) , and (20) . 


GENERAL THROUGH-FLOW THEORY 

Because of the enormous mathematical difficulty in solving the 
preceding set of general equations, the essential feature of the three- 
dimensional flow in turbomachines will be investigated with the follow- 
ing two simplifying conditions: (l) The blade rows will be assumed to 

be placed so far apart that the relative flow through any blade row is 
steady. Under this condition, the equations of continuity, motion, 
and energy in the scalar forms are: 


Continuity equation, from equation (l), 

1 S(pW r r) S(pW u ) S(pW z ) 
r Sr + r Sep + SS 


(24) 


or, from equation (23), 


± S(h* r_1 e“ As W r r) 

r Sr 


1 

+ — 

r 


S(h 


*r-i e -As 


w u ) S(h 

+ 


JL 
¥T-1 


■As x 
e WJ 


Sz 


= 0 (24a) 


The three equations of motion, from equations (14), 


W, 


_u 

r 


"S(V u r) 

Sw r " 

+ w z 

/aw r 


[—sr- 

“ Sep _ 

Isr ■ 

■ w) 


W. 


_r 

r 


S(V u r) 

¥ 
L 

- w z 

(x _ 


[““Sr 

Scp”_ 

\r Sep 

Sz / 


Si 

' SS 

l Si 


T Ss 


r Sq> + r Zq> 


- W T 


Sw r Sw z 

ss - ■ sr. 


1 Sw z 

+ W u li 


Sw. 


u 


r Sep Sz , 


Si 


Ss 


SS + T SS 


and the energy equation, from equation (20), 


Q 


= t( 


w. 


w u Se 


Ss 


r SF + TS^ + Wz SS 


Si 


Wu Si 

= Wr Sr + “r"Scp +Wz SS 


Si 


(25) 


(26) 


(27) 


(28) 
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(2) The Hade is so thin and the number of blades is so large or 
the pitch is so small compared with the blade chord that the terms con- 
taining the circumferential variation of the velocity components in 
equations (24) to (28) are much smaller than other terms in the same 
equations and can therefore be neglected. The circumferential vari- 
ation of pressure or enthalpy is, however, preserved in equation (26) 
by the introduction of a blade force F, which may be considered as 
either due to the pressure exerted by the blade surface on the cir- 
cumferentially thin gas stream or as a circumferentially averaged blade 
force on the gas stream between two blades. This simplification is 
first introduced by Lorenz (reference 8) in order to follow the flow 
along a given surface. The physical concept involved is clarified by 
Stodola (reference 9). Ruden (reference 10) further proves, for incom- 
pressible flow, that the solution so obtained will give an average 
value in the circumferential direction for a finite number of blades, 
provided that the departure from the average value is small. For turbo- 
machines having relatively thin blades of moderately high solidity, 
this solution can indeed be taken as that for a relative stream sur- 
face, which is about midway (based on mass flow) between two blades. 
Because the circumferential variation of pressure and density is con- 
sidered in equation (26), it is better to refer, to the present theory 
as a "through-flow" theory (following Ruden), or "large-number- of -thin- 
blades" theory, instead of "axially symmetric" or "inf init e-number- of - 
blades" theory. (Differentiating velocity components obtained in this 
theory and combining them according to equation (16) does not give the 
true vorticity.) 


or 


With this second condition, equations (24) to (28) become 

1 c>(pW r r) d(pW z ) _ 

r 5r 5z 


(29) 


1 d l ,*r-l -As* ] d (,*r-l -As* . n 

\h e W r ry + ^ \h e W z / = 0 




(29a) 


W, 



di 

d(V u r 

577 

= " 5r + "7 

dr 

dH . 

V u ^(V u r) 

m 5s 

5? + 

r dr + 

T S?- 


+ T 5? + F r 


+ F, 


(30) 
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and 


or 


D(V„r) S(V„r) 


W T 


' U ' 

Dt 

- w r 

dW r 

**z) 

sr - 

5F7 


Sh 

~ — 

^ + 


u * ' d(V r) 

+ W, = F„r 


' z "5z ~ *'u* 


Si . w u d(v u r ) m 




+ T 5S + Fz 


Sh , \ S (V) . m Ss 
r oz 


+ T + F , 


Q = T 5? “ T ( w r sf + W z 55 


Q = 2i — 4 - y 

Dt r cir z 3Z 


(31) 


(32) 


(33) 


(33a) 


The circumferential enthalpy or pressure gradient in equation (3l) 
is replaced by the term F u . The corresponding effects in the other two 
directions are represented by F r and F z in equations (30) and (32), 
respectively. The vector F is perpendicular to W or related to the 
shape of the blade or flow surface by the following equations: 


F r W r + F U W U + F Z W Z — 0 


(34) 


or 

F r dr + F u r dcp + F z dz = 0 


(34a) 


Only six independent equations exist in the preceding equations : one 

continuity relation, three equations_of motion, one energy equation, and 
one orthogonality relation between F and W. For example, equation 
(33a) can be derived from the equations of motion, equation of energy 
(33), and equation (34). (See reference 4.) 


DIRECT AND INVERSE PROBLEMS 

In a direct problem, the blade surface or flow surface is considered 
to be given by the equation 


S(r,0,z) — C 


(35) 
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and two relations among the force components are obtained by relating 
them to the partial derivatives of the function S: 



dS 

pSs 

r Se 


(36a) 




5s 

5? 

rss 

r 50 


(36b) 


Equations (29) to (34) and (36) therefore provide eight independent rela- 
tions that completely define the direct problem involving the eight pri- 
mary variables Vf r , W u , W z , F r ; F u , F z , s, and I or H. (The 
heat transfer Q between the gas and the blade is only important in 
the case of cooled turbine blades and is considered to be given by cool- 
ing considerations.) 


In the inverse problem, the blade surface is to be determined, 
which means that the two relations among the force components as given 
by equations (36) are not available. This unavailability does not mean, 
however, that the designer has freedom to prescribe, the variation of 
two variables or two conditions among the, variables throughout the blade 
region because, in order that the differential equation (34a) will lead 
to an integral blade or flow' surface of the form of equation (.35) , the 
following necessary (and sufficient) condition of integrability must 
be satisfied (reference’ll): 


F-VXF = 0 


(37) 


which, for the present case, reduces to 


5 ftz \ _ S / for \ 

Sr \F u r j “ Sz \F u ry 


(37a) 


This condition of integrability was first pointed out by Bauersfeld as 
early as 1905 (reference 12), but is neglected in many recent investi- 
gations. In effect, it gives a restriction to the velocity variation 
that a designer can specify through the force terms in the motion 
equations (30) to (32) . 


Whereas in the direct problems two conditions are obtainable from 
the given surface, in the inverse problem one condition on the surface 
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must "be satisfied. Hence, the designer can specify, in the inverse 
problem, only one relation throughout the blade region, which can be a 
variation of one thermodynamic quantity of gas (V u , V z , p, or h), one 
relation on the blade surface (untwisted or radial blade element), or 
one relation among the gas properties (constant Mach number), and so 
forth. The variation prescribed should, of course, be reasonable so 
that the solution exists. 

After this one relation is prescribed by the designer, the solu- 
tion of the inverse problem is quite similar to the direct problem. 
Among all the equations to be satisfied, equations (29) to (34) are 
common to both problems. In addition to these six equations, equa- 
tions (36) are available in the direct problem and equation (37a) and 
the one relation prescribed by the designer to be satisfied are avail- 
able in inverse problem. If, in the design, a condition on the blade 
surface is specified such as the blade design in which all blade ele- 
ments are radial, F r is prescribed as zero by the designer, and the 
integrability condition equation (37a) leads to a simple relation be- 
tween F u and F z . These same two relations are also directly given 
by equations (36a) and (36b), respectively, in the direct problem, and 
consequently the solution for the two problems is exactly the same. 

In other cases, however, the solution of the two problems is a little 
different . 


PRINCIPAL EQUATION FOR TWO MAIN GROUPS OF DESIGN 

In this and the following sections, a general method of solution 
for both the direct and the inverse problems will be described. From 
the preceding equations the through- flow considered herein is essen- 
tially described by the equation of continuity (29) and the equation 
of motion either in the radial direction (equation (30)), or in the 
axial direction (equation (32)). Except in the case of low-speed 
centrifugal impellers, it is always advantageous to use equation (30) , 
because F r is either zero in high-speed centrifugal or mixed-flow 
impellers or relatively small in axial-flow bladings. Also, either 
cffl/c)r or bl/br is usually equal to zero, and B(V u r ) fbr is often 
given. This choice is used hereinafter. If it is desirable to use 
equation (32), equations can be developed in a similar manner. . 

The combination of equations (29) and (30) is carried out by the 
use of a stream function, which is defined as follows and satisfies the 
continuity equation (29): 
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> 

) 

1 

1 


% 


Sijr 

5? 


rp*W z 


rtf 


1 

,* 7-1 


As* 

e W„ 


(38a) 


dijr 

dz 


-rp*W r 


-rh* r-1 e _As W T 


From equation (37), 


1 

rh *T-l e -As* 


sw z 


1 _ _1_ dh + ds*\ St 
r a 2 Sr dr y Sr 


-rh* r ‘ 1 e- Aa * 


3w r _ a 2 t 

sr = s? + 


1 Sh Bs*\ St 
12 Si + Si _ /SS 





(38b) 


(39) 


The succeeding development is a little different for the two main 
groups of designs to be considered. In the first group, the variation 
of the angular momentum of the gas about the axis of rotation is pre- 
scribed by the designer; that is. 


V u r = fp(r,z) 


(40) 


is given. Among this group are the free- vortex type in which fj_ is 
just a function of z, the more general solid-body-rotation type, the 
symmetrical- velocity-diagram type, and others. In the second group of 
designs, the following relation between tangential and axial velocity 
is prescribed by the designer: 


w u 

p = tan p = — = f 2 (r,z) (4l) 

w z 

Among this group are the common blade design for higtf-speed centrifugal - 
and mixed- flow impellers in which all blade elements are radial with 
p = r fj(z), the less general design with p = tan p = f^(z), which 
gives a practically untwisted blade and is most suitable for cooled tur- 
bine rotor, and others. The principal equation will now be obtained in 
a form most convenient for these two main groups of designs. 
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Group of Designs in which Equation (40) is Specified 
From equations (9) and (38), 


, T w u . arr* 1 -2,*r-l 2As* 

h = 1 - — + — - 2 r h 6 


s*) 2 + M 2 


Differentiating with respect to r and z yields 


Sh 

S7 


a 


a 2 -W m 2 


Sr Wu Sr 


Si rT u "u . 2 


- W m 2 (- 


1 . Ss* 

x + zr. 


r-V S5| 

Sr 2 


z — _ " r SF^S 


Sh. _ 

a 2 -W r 2 


a 1 - 


m 


Si y ^ W u w 2 Ss* 

Sz Wu SF" Wm ST 


r» z - w r gi 


Substituting equation (43) into equation (39) and adding yield 
1 

1 -As* 

* V 


, 2 TT 2x ,*7-1 -As*/^ W z sw r\ ,2 „ 2x S 2 * „ S 2 * 

(a " w m )rh e I gj- - ^ J - (a -W r ) -3 - 2 W r W z + 


, 2 2^ S 2 i]f ( a 2 Si tt 5W u 2,2 Ss*\ Sf 

(a -W z 2 ) + - “ - 57 + % gj- - ® r + a g- J gj + 


Sz^ 


SW- 


u 


si + Wu sf + 


Ss* \ St 
Sz / Sz 


= 0 


Substituting equation (44) into equation (30) and dividing by o. c 
the following principal equation: 


(42) 


(43) 


(44) 

yield 
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v /\ a i 2 * 


' 2 tr snfe + ^ - -jr) aS + ^ s? + ”1 5J 


M r W 2 a 2 » + . 5*g_) + L, ^t+N-, 


0 (45) 


•where 


Li = 


1 _1_ 
r + a2 


Si tt SW u 2.2 Ss~* , 
^ + W u ^-cor + a S7 - + 


a 2 -W m 2 f SI + S +F +T |s 

Sr r Sr r or 


W 7 


Ni 


i_ ( Si 
a 2 \" Sz 


+ w. 


Sw n 


u 


u ST* 


+ a 


2 Ss* 

sr 



With the variation of V u or W u prescribed by the designer, the 
meridional velocity components are determined by the principal equa- 
tion (45) . The other equations are used to determine various terms 
involved in the coefficients Lp and Up. From the coefficients of 
the second derivatives, the pr incipal equ ation (45) is hyperbolic when 

the meridional velocity W m = ^ W r 2 + W z 2 is greater than the speed of 

sound, and elliptic when the meridional velocity is less than the speed 
of sound. For the hyperbolic case, the method of characteristics can 
be used and will be discussed later. For the elliptic case, it is con- 
venient to put the principal equation in a slightly different form, as 
follows: From equation (38), 


i <*£ d(lLoge p*) * SW r _ S(log e p*) ^ ' 

Sr ^ r 2 r Sr Sr Sr “ rp Sz ^ z 2 Sz Sz 

(46) 

Substituting equation (46) into equation (30) results in 
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In this form, all terms except the first three are taken as constant 
during the successive improvement of \[r values throughout the whole 
region in the numerical solution so that the coefficients in the differ- 
ence equations of ^Jr T s will not vary during successive cycles. 

Equation (45) or (47) is then the principal equation for this group 
of designs to be used for a meridional velocity greater or less than 
the velocity of sound, respectively. The process involved in solving 
this principal equation, together with other equatidns in the inverse 
and direct problems, is as follows: 

In the inverse problem, V u or W u is given by equation (40) . 
Equation (3l) is first used to compute F u . The energy equation (33) 
is then used to determine the variation of s along the streamline. 

The variation of I along the streamline is obtained from equa- 
tion (33a). Equation (32) is used to compute F z and F r is then 
obtained by integrating equation (36a) along a constant r line: 



where F r = 0 at z = z Q . The solution is then carried downstream by 
equation (45) in the hyperbolic case, whereas successive sets of 
improved values of t are obtained throughout the region in the ellip- 
tic case. The quantities W r and W z are then computed from equa- 
tion (38) . 

In the direct problem, equations (36a) and (36b) are given. It is 
most convenient to obtain W u from equation (34) as follows: 

W » - - W r + »*) < 49 > 

The quantity F u is then computed from equation (3l) and F z and F r 
are obtained from equations (36) . Equation (33) is used to determine 
the variation of s along the streamline and equation (32) to deter- 
mine the variation of I or H. Equation (30) is then used to solve 
for \|r as before. 
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Group of Designs in which Equation (4l) is Specified 


For this group of designs, it is necessary to combine W u into 
W z according to equation (4l) as follows: Substituting equation (4l) 

into equation (30) yields 


a^ 2 ) !r - lr + + IK + 20,1 + r z {- 1 * 1 1 + - 0 (50) 

Instead of equation (42) , 


h = I + 


oo 2 r 2 


|r- 2 h* ^e 2 * 8 


(l + u 2 )(^ 


5? 


2 + 

+ \5Z. 


(51) 


Differentiating with respect to r and z, combining with equa- 
tion (39), and substituting into equation (50) give the following form 
of the principal equation: 




J 2 3? + N 2 


where 


1 . 5s* _1 fbl 

a/- 


' v + 5F" " 12 IS? + ® r ' w z - ^ 


L 2 = (l+p 2 ) 

a 2 -W 2 f\i , Sp^ , a 2 -W 2 ( 5l , „ 5s , „ , 0 T r \ 

— Kf + sf; + ^ l- 57 + T a? + F r + ^uj 


(52) 


w Ss* i Si w z Sp 
W 2 = ^'-^5I + ^-^^ 


Equation (52) becomes hyperbolic when the relative velocity is super- 
sonic, and elliptic when the relative velocity is subsonic. For sub- 
sonic velocity, a more convenient form of this equation for computation 
is obtained by substituting equation (46) into (49): 


vi 
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St dpi 

+ 


(53) 


The process involved in solving the principal equation simultane- 
ously with other equations in the direct and inverse problems is essen- 
tially the same as that given for the previous group of designs. The 
only difference is that W u is now computed from equation (4l) and, in 
the elliptic case, is to he reevaluated after successive improvements 
of ijr. In the design where radial "blade elements are employed, the com- 
putation is considerably shortened with F r equal to zero, which shows 
the advantage of using the cylindrical coordinate system for this prob- 
lem. (in the customary method employing a coordinate system along the 
streamline and normal to it, the blade force along the normal is not 
zero.) In this design, the process for the direct and the inverse 
problem is exactly the same. 

The use of p = f3(z) in the design will lead to a blade close to 
the untwisted type if the hub -tip- radius ratio is not too small. In 
the direct problem with such a blade given, this relation may be used 
or the flow may be more accurately obtained by using the equations 
given for the first group. 

The different character of the principal equation for the two 
groups of designs considered is interesting. The character of the prin- 
cipal equation depends on the variation prescribed in the design or 
considered as given in the direct problem. This fact may be utilized 
to solve some flow problems in which the flow in some region is slight- 
ly supersonic. When a tangential-velocity variation is given, the 
equation for the whole region may be elliptic. 


GENERAL METHOD OF SOLVING PRINCIPAL EQUATION 
Elliptic Case 

For convenience of discussion, the fundamental equation (47) or 
(53) of the two groups of designs can be written in the following gen- 
eral fcnn: 
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Si|r 

5? 


S 2 t 

Sr 2 


K d* dfj 

r Sr Sz 2 


?(■ 


St sp St sp\ 4 

J S? S7 + Si SJj + L p 


r 2 p* 2 M = 0 (54) 


The equation is nonlinear even for incompressible flow. The equa- 
tion may he more conveniently rewritten in a linear form as 


S 2 t K St S 2 t _ w 
J ^2 - F 3? + - H 


(55) 


where 




if St Sp St Sp 


2 *2 

L P* + ^** 

s ? 


(56) 


and is evaluated from an approximate solution at the start of the cal- 
culation and from the t and p values obtained in the previous cycle 
during the calculation. For simple boundary shapes and simple func- 
tions of J and K, it may be possible to find a Green* s function 
G(y>z,r|,£) with its proper properties so that the solution of the prob- 
lem can be written in the following form: 

*(r,z) = /fair, z,n,C) N(n,0 ^ (57) 

For example, for flow with the design in which the tangential velocity 
is prescribed, the principal equation takes the form 


S 2 t 

Sr 2 


Sijr 

Sr 


S 2 t _ 
Sz 2 


N 


(58) 


If the boundary walls are cylindrical surfaces, the total- enthalpy and 
entropy are uniform, and the tangential and the radial velocity both 
vanish at the inlet and the exit, the radial variation of i|r at the 
inlet and the exit are the same and can be subtracted from ^Jr, which 
results in 


d 2 (t-ti) l d(t-*i) S 2 (t-i 1 ) 

Sr 2 " r + dz 2 




( 59 ) 
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which can he written as 

+ a , aL (tlE = I 

Sr 2 V r 7 Sr ^ r 2 y Sz 2 V r / r 


(60) 


The quantity 




is zero on the boundary and the 


Green's function is available from a similar equation 
When the Green's function of reference 3 is used, the 
is 


corresponding 

in reference 3. 
solution for \|r 




(61) 


If G is tabulated at several values of r, ijr-ti can conveniently 
obtained by a numerical double integration process on a punch-card 
machine. For curved boundary walls in the meridional plane, it is 
necessary, in this method, to use the technique of conformal trans- 
formation to render the given boundary shape into a rectangular one . 
Inasmuch as this process involves a numerical solution of the Laplace 
equation with the given boundary shape, it is found better to solve 
directly the given equation (55) with the given shape by the numerical 
method. Furthermore, this solution will be the only choice in the 
general case where J and K are not equal to 1 or the boundary con- 
dition is more general, which makes the task of obtaining the proper 
Green's function a very difficult one, if not impossible. 


In order to solve directly the given equation (55) , a general 
numerical differentiation formula for first and second derivatives with 
function value given at unequally spaced grid points using second- and 
higher-order polynomial representation is required to give conveniently 
and accurately the finite-difference expressions at the grid point near 
the curved boundary, which is done as follows: When the value of any 

quantity q is known corresponding to a number of unequally spaced 
values of the independent variable x, denoted by Xq, x, . . . x n , 
the variation of q with respect to x is most conveniently expressed 
by a Lagrangian polynomial of the degree: 


q(x) 


n 


E nn+1 

j=0 


(x) 


x-x4 H n+1 (x4 ) 


o. 


n 


R 


(62) 
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where 

rc ii+l( x ) = (x-XqJCx-X!) . . . (x-Xq) 

n n+l( x j) = ( x j- x o)( x r x l ) ' • * • 

and 



n n +i(x) 
(n+l) ! 


q.( n+1 ) (u) 


(63) 

• ( x j _x n) ( 64 ) 
(65) 


where V lies "between xq and x n . The successive derivatives of q. 
with respect to x at any point x can be expressed as (reference 13) 


n 


i^-z: 

5=0 


qp n n+ i(x) 1 

n n+l( x j) X -Xj X_Xk 

k/j 


+ i 8 


= 21 fz 


3= 0 


3 n n+ 1 (x) A/ i i 


rc n+l( x j) x-xj ^ x-x k f- 


■< x-o 


k+1 “ X * 

k^j \ l^j 


+ n R 


and so forth , with 

& - ^Tjr i (n+1) M + 2 ^ +2) (”) %*■■■• 


T^iyr 


dx 


(66) 


(67) 


( 68 ) 

(69) 


The summation operation is very easily performed when x is a grid 
point, because most of the products vanish. At these points, it is 
convenient to write 


n 




(70) 


The differentiation coefficients B and the coefficients of the deriv- 
ative in the first or second remainder term have been explicitly 
expressed in reference 13 in terms of the spacings between the succes- 
sive grid points for general nonuniform spacing throughout and for the 
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special case near a tapered "or curved "boundary where only the first 
or last spacing is different from the others, using polynomials of 
second, third, and fourth order. For the special case, these coef- 
ficients have also heen computed for difference ratios from 0.1 to 
1.29 in intervals of 0.01 of the distance between the boundary and 
the nearest point and the other spacings and are given in refer- 
ence 13. For spacing lying between tabulated intervals, the inter- 
polation coefficients given in reference 14 can be used to obtain 
the required values of B. 

In the present fluid-flow problems, it is necessary to cover a 
large region in order to reach the boundary condition that is always 
given at stations far upstream and downstream of the blade row. In 
order to reduce the labor of computation, it is desirable to determine 
if the number of grid points required for a given accuracy can be 
reduced by using an order of polynomial higher than the customary 
second order. A study of the remainder terms (reference 13) and actual 
experience in the present problems show that, in most cases, the use 
of fourth-order polynomial will reduce the necessary number of grid 
points to less than one-quarter of that required by the second-order 
polynomial. In setting up the grid pattern, it is always desirable 
to map the flow region in such a manner that the distance between the 
boundary and a point next to it is not too small compared with the 
other distances, because the differentiation coefficient becomes very 
sensitive to the small ratio. If the small ratio cannot be avoided, 
it is best not to include these points in the calculation. 

When the grid pattern and the order of polynomial representation 
have been selected, the coefficients B at each point can be obtained 
from the table given in reference 13. Then the differential equa- 
tion (55) at any grid point whose i|r value is \|4 (fig- l) is replaced 
by the following algebraic equation: 


j=0 J r 1 k=0 


(71) 


where ^ and ^ denote the values of i| r along constant z and 
constant r lines, respectively. 


The values of i|r along the huh and the casing walls can be arbi- 
trarily chosen, with the difference proportional to the mass flow 
between them. At the first station on the left 1-1 and at the last 
station 2-2 on the right (fig. 2), however, the, \|r values are unknown. 
The boundary condition at stations 1-1 and 2-2 is, usually, that the 


/ cnz 
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flow is parallel to the hounding walls. When the huh and casing walls 
are horizontal at the inlet and exit of the machine and the stations 
chosen are far enough out, the value tyb required at a point S z 
distance away from the i|r a point is also equal to Whether or not 

the inlet and the exit stations are chosen far enough out will *he 
indicated hy the variation of \|f's along these stations obtained in 
the solution. 

After equation (71) is obtained at every interior grid point, a 
number of methods can be used to solve the set of simultaneous alge- 
braic equations. For hand computation, the relaxation method developed 
by Southwell has proved to be superior to others for this type of 
equation (references 12 and 15 to 18) . If the fourth-order polynomial 
representation is decided, the calculation can be carried out in two 
steps by using only the five main coefficients in early stages and 
then the residuals are recomputed and relaxed by all nine coefficients 
(reference 13) . 

If a high-speed, large-scale, digital computing machine is avail- 
able, the set of equations can be solved either by an indirect or 
direct method. In the indirect method, the straight iterative method of 
Liebmann is used, wherein ty 1 is solved at each point from its sur- 
rounding values according to equation (7l) and the process is repeated 
until the change at any point is no longer significant. This method 
is simplest to set up, but is slowest. A better method is to set the 
machine to perform a simple relaxation process by computing the resid- 
ual at each point and relaxing according to a fixed relation with 
respect . to the amount of residual just found. 

A much quicker machine method, especially when a number of solu- 
tions with, for example, different inlet Mach numbers are to be found 
with a given geometrical shape of the problem, is the direct method 
that solves the set of simultaneous equations by a matrix process. 

The details of a matrix method > which fully utilizes the great number 
of zero elements of the original coefficient matrix, is given in refer- 
ence 13. 


Hyperbolic Case 

In the hyperbolic case, the problem consists in solving i|r from 
the following principal equation, written in a general form for the 
two groups of designs: 



28 


KACA TN 2302 




£ 

ciz 2 


(72) 


with the initial condition that \|r and its normal derivative are given 
on a curve that is not a characteristic line. 


From equation (72), the equation of the characteristic line is 



The slopes of the characteristic lines in the r,z-plane are 



(73) 

(74a) 

(74b) 


When the point on the r,z-plane moves along the characteristic 
curve l, corresponding to a small change dz in z, the change in 
r is dr = dz. Because of these two small changes, the change of 
any quantity q is 


= § dz = S? dz + S? x i dz 


(75) 


or 


dq _ Sq , Sq 
dz Sz 1 Sr 


Hence, along Z-^, 

A (V) + a (k) = 

dz \3zy dz \ozJ 1 or \oz J 


s 2 t + x 

^ z 2 1 dz dr 


d / djr 
dz \ dr, 


(1 j) + X 1 SF 


S 2 \jr . 

dz Sr 


d 2 fr 
1 Sr 2 


(76) 


(77) 

(78) 


* 


; on? 
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From equatio ,) and (78), 

d 2 t _d_ _ x 

dz dr dz \Sr / 1 ^ r 2 

d 2 t d_ /dijr) x , x 2 d 2 T]r 

^ z 2 “ dz \dz/ " 1 dz v5r/ 1 ^ r 2 

Substituting equations (79) and (80) into equation (72) yields 

T d /d 


(79) 

(80) 


dz \dz/ + ( 2K “ jX l) ^ (H) + (V- 2ia L + L)^| + M^| + N^|=0 (81) 
By virtue of equations (74a) and (74b), equation (8l) becomes 


Similarly/ along the second characteristic line Zg, 

f m + f (§*) + a It + 5 |t , o 

dz \dz / 1 dz \dr/ J dz J or 


(83) 


Starting from two points a and h a short distance apart on the curve 
where the initial condition is given, equations (74a) and (74b) give 
the tangent to the characteristic curves at these two points and equa- 
tions (82) and (83) give the new value of cty/Sz and S^/5r at the 
point of intersection c of the two tangent lines. The auxiliary 
equations corresponding to the particular problem are then used to 
determine other pertinent quantities at the point c. This process is 
to be carried step by step downstream. The method is the same as for 
ordinary two-dimensional rotational flow. (For details of calculation, 
see reference 19 . ) 


APPLICATION TO TYPICAL DESIGNS 

The following sections include a brief discussion of the manner 
in which the fundamental and auxiliary equations reduce to particular 
forms for several typical designs. In actual computations, it is 
always desirable to render all quantities dimensionless. A convenient 
system is to divide r or z, W or V, p, T, s, H, or I, 
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F, *, and 03 by r t , U t , P T ^, U t 2 /k, R, U t 2 , U t 2 /r t , r t 2 U t , 
and U^/r-j., respectively. These dimensionless values are used in 
the following equations. 


Free-Vortex Design 

In the free-vortex design, the variation of V u r is prescribed as 
a function of z only. With a free-vortex flow, the total enthalpy at 
any point z is simply related to the inlet value by 


H*-H? = 


J j _Q_ 
i u t 2 


dt + o?[(V^) r (V^)J 


(84) 


where (V* r *)i i- s a constant. If is uniform with respect to r 

and SQ, /dr is zero, dH/dr will be zero everywhere; but dH/dz is not 
equal to zero in the rotor, whereas dl/dz is zero and it is there- 
fore convenient to use the system of equations involving I. The 
principal equation is then 


d 2 ** ± 

d r *2 ' r * dr* dz*2 


1_ (b£ d^ df * dp^\ r* 2 p* 2 

p*Vdr*df* dz*dz*J 

dr* 


•* 2 P* 2 ( F * + m* $2* \ _ 0 

d£ ( F r + T to *)- 0 


(85) 


The auxiliary equations are 


*: * 1 dty* ^^ u r | 

dr* dr* 


F u r 


( 86 ) 


* _ Wu_ ^( v £r*) 
F z " r * 


dz* 


V. * 

rp* 2® _ 

vLr 


d*’ 


dz* r * 2 p* 2 dz* 


dV 


1 d\jf* d 2 \[r* 


i /d r d P * d»* d P * 


dr* 2 r* dr* dz* 2 P* \dr* dr* 


dz* dz*> 


(87) 


F * = 4 (if?) dE 


‘( 88 ) 
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?* 


* 



The computation can he started with an assumed value of i|r con- 
forming to the boundary shapes. In the early stage, it is advantageous 
to omit equations (87) and (88) and to use the following approximate 
formula of F r , which is obtained by solving the two equations by 
assuming ch|r/cir constant and neglecting small terms: ' 


^ 2 ^Ofr*) 

F„ = 




II 

Jzq 


r* d£ 


(90) 


Density can be obtained from the \|r derivatives as follows: First, 

2 

equation (89) is written as a = [l - cp/cr) , where 


a = p^e 2 * 3 


*2 


£_ c»* 2 r* 2 W u 

H* + 2H| 


2H 1 


2 

r-i 


cp = 



+ f^!f 


U z *y J 


e 2As* (2r* 2 H^) _1 I * ■ co * 2r * 2 Wu 


*2 


y+1 

r-i 


.H* 2H* 


2H* 


Second, either cp is computed for a number of values of a and a curve 
is plotted, or a is obtained for a number of values of cp through 
iteration of the preceding equation and a table obtained for equal 
intervals in cp. In either case, after the t|r derivatives are obtain- 
ed, cp is computed and o or p is then obtained from either the 
curve or the table. 


Design Based on Symmetrical Velocity Diagram at All Radii 

Generalized for compressible flow, the symmetrical velocity dia- 
gram at all radii is defined as follows (reference 4): (See fig. 2 for 

station notation.) 
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«*(*&* )* a(V*r*) e * 


to* 




= 2r> 


(91) 


If it is desirable to maintain a constant total state along the radius, 
the variation of V u r with z is such that the same amount of work is 
done along all streamlines: 


(V u r ), - (V u r ) b = f(zj) 


Hence, 


^(Vu r *)j d(<r*) b 


Combining equations (91) and (93) yields 

w***h. r * 

~W~ b 


= 0 


(92) 


(93) 


(94) 


Hence, 


*2 *\ , T7 * 
r b (V )e,t-( V u r )b,t 
(V )b - -2 i T - " 

,„** x r b (V )e,f( V u r )b,t 

(V u r ) e = — " + 2 “7 


( V J r *)j = ( V u r *)b + [( V i r *)j,t - ( V I r *)b,t] 


(95) 
■ (96) 

(97) 


In this design, the quantity Sl/^r is not zero in the rotor, 
but with constant work input, bft/br is zero. Using the system of 
equations involving H is therefore convenient: 

d 2 t* 1 5 2 t* 1 / dfr* 5p* 

5r* 2 ’ r* Sr* Sz* 2 “ p* Ur* Sr* Sz* Sz* 


*2 *2 
r P 

S\|f* 

Sr* 


-fc 1 
F u r = -?-* 


V* S ( vjr * ) * 


* Ss* 


S^~ + F r + T S? 


= 0 


a/aofr*) 5 t »s(v; r *) 

Sz* Sr* + Sr* Sz* 


(98) 


( 99 ) 
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*?- 


T ft -v / % \ 

W u S ( V u r ) 


dz* 


- T 


* Ss 

s ? 


- (JO 


•s * 

* Bz* 

a £ 

Sr* 


S(V*r*) 


Sr* 


*2 *2 
r p 


S£ 

Sz* 


5 2 V* _ _l_ S->[f* S^^* l /Si|r* Sp* Si|f* Sp' 

Sr* 2 r* Sr* Sz* 2 " p* Ur* Sr* Sz* Sz* 


( 100 ) 


* *■ * 

F-v. = F r 
x r u 


,* dr* 

*0 


F n 


F*r* 


<l£ 


( 101 ) 


p* 2 = 


H_ 

H* 


(<r *) 2 


2r* 2 H* 


* 2 

(M_\ 

\Sr*/ 


+ 


/chj^ 


^2 *2* 

2r p 


2 

r-i 


-2As 


( 102 ) 


With only the additional complication that the value of V u in each 
cycle is determined from a knowledge of the streamline in the previous 
cycle , the solution of this problem is obtained in the same manner as 
in the previous case. For a multistage compressor, it is important to 
account for the effect of loss on density rise by including the e”^* 
factor in equation (102) . The increase of entropy can be estimated by 
a knowledge of the polytropic efficiency (reference 4) . 


A nonvortexrtype velocity diagram quite similar to the pre- 
ceding one is the solid-body rotation design that has a tangential 
velocity varying linearly with radius in front of the rotor; that is, 


/ Jk 3k s n 3^2 
(Y u r ) b = Cr^ 


S(V*r*) b _ 




= 2Cr-v- 


(103) 

(104) 


If total enthalpy Is to he constant along the radius through the rotor, 
equations (92) and (93) also apply. Hence, 

(Vur*)j = Cr £ 2 + [(vjr*)j,t - (vjr*)h,t] 


(105) 
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and 

(vjr*) e - Crf + [(V*r f ) e;t - (<r*) b/t ] (106) 

Comparison of these equations with those in the previous design is 
interesting. In the previous design the change of whirl through the 
rotor is distributed evenly in the whirl in front of and behind the 
rotor, whereas in the present design it is completely put into the 
whirl behind the rotor. Except for this difference, the calculation is 
quite similar to that given in the previous design. 


Designs Involving Untwisted Rotor Blades 


Untwisted rotor blades are desirable because of simplicity in manu- 
facturing, and seem to be the most practical design for cooled-turbine 
rotor blades. They can be efficiently used if the stator blades are 
designed to fit them. If the blade is not too long, it can be designed 
on the basis that p = W- U /W z is a function of z only. The principal 
equation then takes the following form: 


(1+h 2 ) 


s 2 t* 

Sr * 2 


Sr* 


S 2 t* 

Sz* 2 





st^ s £ 

Sz* Sz* 


2pco*r* p* + 


r* 2 P* 2 l* 

St* V r 

Sr* 


Si* 

Sr* 


+ 


i* ^ 

' Sr* 


= 0 


(107) 


Radial- and Mixed-Flow Impeller with All Radial 


Blade Elements 


The speed of rotation of the rotor can be increased by having all 
blade elements radial. With F r = 0, the integrability equation (37a) 
gives the result that F z /F u r is a function of z only. Thus, when 

-| = r* f 3 (z*) (108) 
3 


equation (41) becomes 



F* 


\x = tan (3 


r* f 3 (z*) 


(109) 
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The principal 


ion is then 
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l+(r* f 


l+(r* f 5 ) 2 d 2 V* 


3 ' J^ r *2 


dr* dz : 


,*2 


p*L 


+ „*2 r 2\d$t df£ dtf^_ dp* 

3 dr* dr* dz* dz* 


- 2vo?r*p* - 


r* 2 p* 2 /dl* * _ 

at* \dr* dr*y 

dr* 

( 110 ) 


With F r = 0, use of the three auxiliary equations to compute F r is 
unnecessary. 

Equation (110) is further simplified with di/dr equal to zero if 
the inlet flow is of the free-vortex type or has no whirl. 


Simple Approximate Correction for Blade Thickness 

If the "blade is not quite thin, it is desirable to add a simple 
approximate correction factor b in the definition of stream functions 
of equation (38) as follows: 

1 

, cMlr * TT tJfT-I -As^tt 

gj = rp ¥ 2 = rl1 e w z 

i 

b Sz ” ~ rp*W r “ “ rh*^“ 1 e“^ s *W r 

A good conception of this thickness correction factor can be 
obtained by analyzing the effect of blade thickness on the specific mass 
flow along the mean streamline (based on mass flow) in two-dimensional 
cascades. Yet unpublished calculations made for a typical subsonic tur- 
bine cascade and two supersonic compressor cascades show that the 
specific mass flow pW z on the mean streamline is about 4 and 10 per- 
cent higher than that given by one -dimensional calculations correspond- 
ing to the same reduction in channel area for the subsonic and super- 
sonic cascades, respectively. The influence of the blade thickness 
also extends a short distance upstream and downstream of the blades. 

The shape. of the mean streamline is also seen to follow .approximately 
the mean channel line, (but with lower curvature) . When this correc- 
tion factor b is used in equations (ill), all the equations previously 
obtained should be modified by replacing p by p/b. 


(111a) 

(mt) 
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CONCLUDING REMARKS 

Equations of motion and energy for unsteady, three-dimensional flow 
of a nonviscous fluid are expressed in terms of gas quantities most con- 
venient for analyzing flow in turbomachines . Entropy change due to 
heat transfer in a cooled turbine and due to strong shock wave in super- 
sonic flow can he taken into calculation. 

The general equations are simplified according to the standard 
assumption for steady through-flow calculation in turbomachines having 
thin blades of high solidity. The problem is completely defined in 
the direct problem with blade shape given; whereas in the inverse or 
design problem, with the inclusion of the integrability condition for 
the blade surface, either one flow variable or one relation among sever- 
al variables can be prescribed by the designer. 

Through the use of the stream function, the continuity equation 
and the equation of motion in the radial direction are combined to form 
a principal equation for the present problem. The principal equation 
contains some terms that are either prescribed or to be determined by 
other equations defining the problem. Two forms of the principal equa- 
tion are obtained for the two main groups of current compressor and 
turbine designs in which either the variation of tangential velocity 
or the variation of the ratio of tangential to axial velocity through- 
out the blade region are given. When the tangential velocity is given, 
the principal equation is elliptic or hyperbolic depending on whether 
the meridional velocity is subsonic or supersonic. When a relation 
between the tangential and axial velocity is given, the principal 
equation becomes hyperbolic when the relative velocity is supersonic. 

A general method of solution for both the elliptic and the hyper- 
bolic cases is outlined. Specific applications of the theory to 
several common types of compressor and turbine employing free-vortex, 
symmetrical- velocity-diagram, solid-rotation, nontwisted-blade, and 
radial -blade-element designs are discussed. A simple correction factor 
for blade-thickness effect is also suggested. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, October 25, 1950. 
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Figure 2. - Boundary conditions and station notation. 



* 


z 


Figure 3. - Characteristic curves 
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